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ABSTRACT 

We have examined in detail the morphologies of seven z ~ 1.5 passively evolving luminous red galaxies 
using high resolution HST NICMOS and ACS imaging data. Almost all of these galaxies appear to be relaxed 
systems, with smooth morphologies at both rest-frame UV and visible wavelengths. Previous results from 
spectral synthesis modeling favor a single burst of star formation more than 1 Gyr before the observed epoch. 
The prevalence of old stellar populations, however, does not correlate exclusively with early-type morphologies 
as it does in the local universe; the light profiles for some of these galaxies appear to be dominated by massive 
exponential disks. This evidence for massive old disks, along with the apparent uniformity of stellar age 
across the disk, suggests formation by a mechanism better described as a form of monolithic collapse than as 
a hierarchical merger. These galaxies could not have undergone a single major merging event since the bulk 
of their stars were formed, more than 1 Gyr earlier. There is at least one case, however, that appears to be 
undergoing a "dry merger", which may be an example of the process that converts these unusual galaxies into 
the familiar spheroids that dominate galaxies comprising old stellar populations at the present epoch. 
Subject headings: galaxies: evolution — galaxies: formation — galaxies: high redshift 



1. INTRODUCTION 

A wealth of observations over the past decade have collec- 
tively converged on the ACDM "concord ance" cosmology , 
including results from high-z supernovae dRiess et al J 1 19981: 
Perlmutteret~aLlll999t iTonrv et alj 120031: IWood-Vasey et all 
2007|), the 2dF Galaxy Redshift Surv ey (iPeacqck et alJl2001h 
Percival et alJl2002t ICofe et al.ll2005h . S PSS (iTeemark et al] 
2004t lEisenstein et all 120051). WMAP dBennett et all [ 2003 1 : 
Spergel et al. 120071) . and weak lensing dHoekstra et alj|20*0*6l : 
Massev et al.ll2007l) . Unfortunately, even given the precision 



of these measurements and the agreement between such a 
wide variety of independent observations, the formation of 
structure on galaxy-size scales remains a challenge for cos- 
mological models. In particular, the formation and evolution 
of the first massive galaxies is still poorly understood. This is 
distressing given that more than half of all luminous matter in 
the local universe is containe d within massive ellip t ical galax- 
ies bu ilt up at earlier times (Fuku gita et alj[l998l: iBell et all 
120031) . 

1 Based on observations made with the NASA/ESA Hubble Space Tele- 
scope, obtained at the Space Telescope Science Institute, which is operated 
by the Association of Universities for Research in Astronomy, Inc., under 
NASA contract NAS 5-26555. These observations are associated with pro- 
gram # GO-10418. 

2 Currently at the University of California Observatories/ Lick Observa- 
tory, Department of Astronomy and Astrophysics, University of California at 
Santa Cruz, 1 156 High St., Santa Cruz, CA 95064. 

3 Also at Cerro Tololo Interamerican Observatory, Casilla 603, La Serena, 
Chile. 



Simple semi-analytic models (SAMs) based on the stan- 
dard cosmological paradigm predict that massive galax- 
ies form rather late in the history of the universe (1 < 
z < 2) from successive mergers of smaller galaxie s (e.g . , 
iThomas & Kauffmannl fl999l: iKauffmann & Haehneld 2000). 
These models, however, have traditionally suffered from an 
overprediction of the number density of both low and high lu- 
minosity galaxies, as well as an inverted color-magnitude re- 
lation (i.e., luminous blue galaxies), and an inability to repro- 
duce the observed color bi-modality in galaxy distributions. 
Feedback from sour c es such as star form a tion and AGN (e.g. , 
iGranato et al.1 120041: ICroton et alj 120061: iBower et alj fZOO® 
seems necessary in order to reconcile these models with ob- 
servations. Results from the latest SAMs predict that, while 
the stars in massive elliptical galaxies formed quickly and 
early on, their assembly into the most mas sive galaxies is 
chara cterized by a late, hierarchical build-up dDe Lucia et alj 
l2006h . 

The predictions of hierarchical assembly can be explored 
through high resolution morphological studies of massive 
galaxies comprising old stellar populations at high redshift. 
Examples of such g alaxies have become increasingly abun- 



dant in the literature ( Liu et al]|2000b 
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phological information has still been relatively sparse. It is 
commonly assumed that spheroids should dominate evolved 
galaxies at high redshift as they do in the local universe. 
However, some recent studies have found evidence for S0- 
type galaxies as well as ex ponential disks among these 
evolv e d galaxies at hi g h z (live et al] 12003b ICimatti et al] 
2004t lYan et al] 12004b iDaddi et al] 12005b iFu et al.1 12005b 
Stockton et al] I2004L 120061 120081) . We present a sample of 
seven z ~ 1.5 "old galaxies" (OGs, by which we refer to 
the age of the dominant stellar population and not necessar- 
ily the age given by the time elapsed since galaxy assembly) 
selected to have both R - K' > 6, and J - K' — 2. We 
have obtained high-resolution morphological data for these 
objects from HST ACS and NICMOS imaging. We discuss 
the range in morphological classes among OGs in our sample 
at z ~ 1.5 and the implications for galaxy formation scenar- 
ios. Throughout this paper we assume a flat cosmology with 
H = 70 km s- 1 Mpc~\ M = 0.3 and fl A = 0.7. 

2. OBSERVATIONS 

In iMcGrath. Stockton. & Canalizol (120071 hereafter Paper 
I), we presented six old galaxy candidates in the fields of 
radio-loud QSOs at z = 1.5, five of which were dominated 
by stellar populations that were already > 1 Gyr old. Here we 
continue our study of these six objects, with high resolution 
imaging from HST. In addition, we include a seventh source, 
J094258+4259.2 ERO, discovered serendipito usly in the field 
of Abell cluster 85 1 by llye et al.1 (2000, 2003). 

Each galaxy was observed with HST at approximately rest- 
frame U and R-bands using the ACS F814W filter (2 or- 
bits) and the NICMOS2 F160W filter (1 orbit). Details of 
the observations are given in Table 1 of Paper 1, except 
for J094258+4659 ERO, which was observed on 2005 April 
18 (UT) with ACS for 5464 s and on 2005 April 19 (UT) 
with NICMOS2 for 2880 s. We used the standard pipeline- 
processed ACS images for analysis. The NICMOS images, 
however, required a more careful processing. We first checked 
to make sure the observations were taken long enough after 
SAA passage (> 30 min.) so that they did not suffer from 
cosmic ray persistence. We then took the calnica processed 
images, which had been flat-fielded and corrected for ampli- 
fier glow, as a starting point and corrected for bias offsets and 
inverse flatfield effects using STSDAS tool pedsky. We gener- 
ated a bad pixel mask from the data quality file, and included 
an additional mask for the coronagraphic occulter which pro- 
duces background ra diation detected in the F160W filter. We 
used L. A. Cosmic (Ivan D okkum 2001) to perform a first cut 
of cosmic rays from the individual images, and then drizzled 
the individual dither positions together to obtain a cosmic- 
ray-free, subsampled image. We performed a series of tests 
with Tiny Tim 4 generated PSFs at the location of our galaxies 
on the NIC2 chip in each of the dither positions to determine 
the best drizzle parameters required to maintain a tight PSF 
as well as recover some resolution in the final image. The pa- 
rameters we adopted were an input drop size of 0.7 convolved 
with a gaussian kernel and an output pixel size of 0.5 times 
the original. 

3. STELLAR POPULATION PROPERTIES 

Stellar populations for six of these galaxies were presented 
in Paper I and were determined from broadband photomet- 
ric fits to S. Chariot & G. Bruzual (2007, private communi- 
cation; hereafter CB07) instantaneous burst models. In that 

4 http://www.stsci.edu/software/tinytim/tinytim.html 



TABLE 1 

Summary of Simple Stellar Populations 



Galaxy 




Age 


Av 


P(x 2 ) 


Stellar Mass 






(Gyr) 


(mag) 


% 


(IO^Mq) 


TXS 0145+386 ER1 


1.4533 


0.90 


0.79 


70 


1.74 


ER2 


1.459* 


2.40 


0.02 


61 


2.02 


J094258+4659 ERO+ 


1.5** 


1.0 


< 1.4 




3.0 


B2 1018+34 ER2 


1.4057 


0.64 


1).94 


6 


1.42 


TXS 1211+334 ER1 


1.598 


1.43 


0.23 


91 


2.90 


AC 15.55 ER2 


1.412 


1.80 


0.02 


89 


1.75 


TXS 1812+412 ER2 


1.290 


1.43 


0.51 


44 


1.56 



* The best-fit photometric redshift is given as no spectroscopy was obtained for 
this source.^ Data are from;Iye et al. ( 2003) and were determined using BC03 
rather than CB07 models.** The weighted best-fit photometric redshift is given 
since no prominent absorption or em ission features were visible in the rest-frame 
visible spectrum of lye et al. 1 2003). 



paper we discussed uncertainties in the models and how the 
age-metallicity degeneracy affects our resulting age determi- 
nations. Extremely deep imaging was obtained at blue wave- 
lengths in order to place constraints on the rest-frame UV por- 
tion of the spectral energy distribution (SED). This spectral 
region is important for age-dating these galaxies, since even 
a small contribution by mass from young stellar populations 
would cause a disproportionate increase in the total light at 
these wavelengths. For this reason, rest-frame UV spectra for 
five of these six sources were also obtai ned and were fit to 
a varie ty of models, includin g the CB07, iBruzual & Charloil 
d2003l) . and lMarastoiil d2005l) models. The results from inde- 
pendent fitting of the spectra and the photometry were consis- 
tent to within ^100 Myr in most cases (see Paper I). Further- 
more, these spectra provide confirmation of the galaxy red- 
shifts, removing a free parameter from the photometric SED 
fitting and resulting in more accurate age determinations. The 
stellar population properties of the remaini ng source, ERO 
J094258+4659.2, were discussed in detail in llye et all (120031) . 
who used both broadband photometry and a low S/N rest- 
frame visible spectrum to place limits on the age of the domi- 
nant stellar population. In all but one galaxy in our sample (B2 
1018+34 ER2), the dominant stellar populations were found 
to be ~1 Gyr or older. A summary of the properties of the 
stellar populations of these sources is given in Table Q] 

The use of instantaneous burst models places a lower limit 
on the age of the dominant stellar population. More realis- 
tic models (e.g., exponentially decaying star formation rates) 
would require star formation to begin at even earlier times in 
order to reproduce the observed spectral energy distributions 
(SEDs), especially the sharp inflection at 4000 A that is the 
signature of old stellar populations. However, given that an 
instantaneous burst of ~ lO n M in stars is not strictly phys- 
ical, it is useful to determine the maximum allowable contri- 
bution to the SED from more recent episodes of star forma- 
tion. For TXS 0145+386 ER1 and B2 1018+34 ER2, where 
we see [O II] emission in the rest-frame near-UV spectra, if 
we assume this emission is due purely to star formation, rather 
than AGN or LINER activity, the rate of current star forma- 
tion is_jjp_2roxm^eI^j^ respectively 
(see lStockton et al.ll2006b IMcGrath et alJl2007h . For galaxies 
without [O II] emission, we can place limits on the amount of 
current star formation directly from the SED modeling. As- 
suming a composite stellar population consisting of an instan- 
taneous burst in which the majority of the stars were formed 
plus a small contribution from an underlying continuous star 
forming population, we can determine the maximum amount 
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TABLE 2 

Composite Stellar Population Properties 



Galaxy 


Age 


Av 


current SFR 


P(X 2 ) 


Stellar Mass 




(Gyr) 


(mag) 


(M yr" 1 ) 


% 


(10 U Af Q ) 


TXS 0145+386 ER1 


1.02 


0.75 


0.85 


47 


1.70 


ER2 


2.75 


0.00 


0.1 


12 


2.15 


B2 1018+34 ER2 


1.90 


0.29 


1.0 


75 


1.88 


TXS 121 1+334 ER1 


2.10 


0.00 


0.2 


52 


3.31 


4C 15.55 ER2 


2.00 


0.00 


0.1 


67 


1.86 


TXS 1812+412 ER2 


1.68 


0.46 


0.1 


47 


1.73 



of current star formation that is plausible given the broad-band 
colors in the SED. We have re-fit each of our galaxies from 
Paper I with this model, allowing both the age of the instan- 
taneous burst and the percent contribution from continuous 
star formation to vary, and we show the results in Table [2] In 
summary, ages for all the galaxies increase (compared with 
the pure instantaneous burst model), and star formation rates 
>lM0yr -1 are excluded for every galaxy except B2 1018+34 
ER2. 

4. MORPHOLOGIES OF OLD GALAXIES AT Z ~ 1.5 

In order to study the morphologies of these seven galax- 
ies, we used our NIC2 imaging data, which have a higher 
S/N than the ACS images and are less likely to be affected by 
any small scale or local variations in the light profile. These 
images should provide the smoothest galaxy profiles which 
will in turn provide the best estimate of the underlying galaxy 
structure. We then compared the NICMOS profiles with those 
obtained from the lower S/N ACS images in order to deter- 
mine if there is any evidence for a difference in the distribu- 
tion of older and younger stellar populations (whose contribu- 
tions dominate the light profile long-ward and short-ward of 
the 4000A-break, respectively). 

We used the program, GALFIT 5 dPeng et alj|2002l) to model 
the galaxy light profiles. These galaxy models were con- 
volved with the Tiny Tim NICMOS PSF, which was driz- 
zled in the same manner as the observations. For all galaxies, 
we determin ed the best fit exponential, r 1 / 4 -law, and Sersic 
dSersidl 19681) profiles. If these simple models did not provide 
an adequate fit to the data, we then tried models with a central 
bulge plus extended disk, or models with a small point-like 
core. 

PSFs were also generated for the ACS imaging in order to 
convolve them with a model galaxy to fit the ACS light pro- 
files. These PSFs required a little more care in making be- 
cause of the distortion in the ACS field. Tiny Tim produces 
PSFs for ACS that are geometrically distorted, such that they 
approximate the PSF at the location of an object in the raw 
data before pipeline processing. In order to use these PSFs 
with the final drizzled and distortion-corrected images, one 
needs to correct these PSFs in the same manner as the data. 
One cannot simply use the pre-distorted PSFs from step two 
of Tiny Tim because the effects of charge diffusion are not in- 
cluded until the final step and the process of correcting for the 
distortion could produce artifacts or other imperfections not 
represented in this ideali zed PSF. We used the IDL wrapper 
for Tiny Tim created by Rhod es et al.l d2007l) 6 in order to pro- 
duce distortion corrected PSFs at the position of our targets, 
which we then drizzled in the same manner as the data to pro- 

5 GALFIT is available at http://www-int.stsci.edu/~cyp/work/galfit/galfit.html 

6 see http://www.astro.caltech.edu/~rjm/acs/PSF/ 



duce the most realistic representation of the PSF for our ACS 
imaging. 

The observed HST ACS and NICMOS images are shown in 
Figure Q] along with the best-fit GALFIT models, the residuals 
from subtracting the model from the NICMOS image, and the 
models without PSF convolution, which are representations 
of the true galaxy light profiles. The observed surface bright- 
ness profiles (Figs. [2] — [H} are plotted at radii increasing by 
intervals of one original (non-subsampled) pixel ((X'075 for 
NICMOS and 0"05 for ACS) outward from the center of the 
galaxy. They are averaged along ellipses of constant position 
angle but varying ellipticity. Model profiles were determined 
using the same set of ellipses as those used for the observed 
data. Adjacent data points and their errors are correlated to 
some extent because of PSF smearing and the fact that the el- 
lipses sample some of the same pixels along the minor axis. 
A summary of the best-fit morphological properties for each 
galaxy is listed in Table [3] 

4.1. TXS 0145+386 ER1 and ER2 

Two EROs with old stellar populations were found in the 
field of TXS 0145+386. From the HST imaging, ER1 ap- 
pears to be a fairly regular elliptical galaxy. Previous re- 
sults from ground-based adaptive-optics (AO) imaging indi- 
cated that this g alaxy may have been a disk galaxy with a 
low Sersic index (Stockton et al. 2006). However, the S/N of 
the AO imaging was not very high, particularly in the outer- 
most regions of the galaxy where contribution from the near- 
IR background becomes important and is both higher and less 
stable from the ground than it is from space. It is at large 
radii where the best-fit exponential and r 1 ^ 4 -law profiles de- 
viate most noticeably from each other. From the higher S/N 
NICMOS data, it is clear that a Sersic profile with an index 
n = 4.29, or very close to a classic de Vaucouleurs ellipti- 
cal dde Vaucouleurs|[r948l) . is the best fit to the observed data 
(Fig. |2]), rather than an exponential disk. The residuals from 
th e fit are within the background noise. 

Stoc kton et alj noted from their AO imaging that there ap- 
peared to be some low surface brightness emission in the out- 
ermost regions of the galaxy, ER1. This is resolved in our 
ACS and NICMOS images and appears to be a faint, com- 
pact, neighboring galaxy. There does not appear to be any 
visible interaction between this faint galaxy approximately 2" 
(17 kpc) to the east of ER1 and ER1 itself, although the two 
galaxies have nearly identical colors. At this redshift, how- 
ever, particularly given the low luminosity of the "compan- 
ion" galaxy, it would be difficult to see the signatures of inter- 
action if the two were undergoing a "dry merger" since these 
features would have a low surface brightness and would be 
dominated by old stars, rather than bright knots of star forma- 
tion. We therefore cannot say whether a merger between these 
two sources is occurring, however, we note that the Hp 160 w- 
band flux ratios of the galaxies implies that any merger would 
be minor and would not significantly disrupt the structure of 
the primary galaxy. 

ER2, the second old galaxy found in the field of 
TXS 0145+386, was unresolved from ground-based, non-AO 
imaging. From the HST imaging it is clear that ER2 is in fact 
an interacting galaxy pair. One of the galaxies in the pair ap- 
pears to be much redder and more extended than the other, 
more compact galaxy. Modeling of the galaxy pair in ER2 
is a bit more uncertain than modeling a single galaxy due to 
the increased degrees of freedom. Nevertheless, the compact 
galaxy in this pair is well modeled by a Sersic profile with 
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FIG. 1 . — HST imaging of old galaxies. Images are 4" on a side, North is up, and East is to the left. From left to right we show the ACS F814W, NIC2 F160W, 
best fit GALFIT model for the NICMOS image, residuals obtained from subtraction of the model from the NICMOS image, and the GALFIT model without PSF 
convolution. 
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FIG. 2.— Surface brightness profile for TXS 0145+386 ER1. Data points in this and the following six plots (Figures 2-8) are from the NICMOS F160W 
and the ACS F814W observations, and are shown with 3<r errors. Key for the models is given in the figure. The bottom two panels show the deviation of the 
data points and other models from the best fit model for both the NICMOS (top) and ACS (bottom) observations. Residuals after model subtraction from the 
NICMOS image are shown to the right of the plots and are smoothed by a Gaussian with width, a = FWHM of the PSF to eliminate pixel-to-pixel variations. 
These residual images are outlined in color according to the appropriate model, as given in the figure key. For TXS 0145+386 ER1, the Sersic and de Vaucouleurs 
models overlap for all but the largest of radii, and provide equally good fits to the data. 



a high index (n > 4) and a small effective radius. We con- 
strained our fit to the n = 4 case, and only allowed GAL- 
FIT to model the second, more extended galaxy. This second 
galaxy is well represented by a profile that has a Sersic in- 
dex intermediate between a classical elliptical and a classical 
disk, although closer to a disk. The best-fit Sersic index gives 
n = 1.86 (Fig.0. The ACS image for this galaxy shows a 
highly irregular structure, and can only be modeled with a rel- 
atively broad and flat profile, similar to that of a low surface 
brightness disk. 

4.2. J094258+4659 ERO 

J094258+4659 ERO is visibly elongated in both the ACS 
and NICMOS images, and appears to be either an edge-on 
disk or SO galaxy with a small central bulge. Comparison 
of the NICMOS and ACS images shows that the stars in the 
bulge appear to be somewhat redder than stars in the disk, as 
expected for typical disk or SO galaxies. It is likely that the 



bulge component formed first, with the disk completing its 
star formation only slightly later. The dominant stellar popu- 
lation, however, appears to be > 1 Gyr old (lye et al. 2003), 
so star formation in both the disk and bulge components must 
have occurred extremely rapidly in order to produce such an 
e volved population by z = 1.5. 

live et all (120031) noted from their Subaru imaging that 
J094258+4659 ERO appeared to be an SO-type galaxy with 
approximately 30% of the total light being contributed from 
a central bulge component. The NICMOS imaging confirms 
this result and the best fit GALFIT model is an exponential 
disk with a 32% bulge component (Fig.|4]i. The ACS image, 
while noisier in the outer regions, also appears to favor a disk- 
dominated galaxy with 25% of the total light contributed by a 
central bulge. 

4.3. B2 1018+34 ER2 
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FIG. 3. — Surface brightness profile for TXS 0145+386 ER2. The profile is centered on the more extended galaxy in this merging pair, and the bump at ~0'/7 
is due to the second, more compact galaxy. Errors are determined by the standard deviation from the azimuthal average, and are therefore larger near the Of.' 7 
bump since this is an asymmetric contribution to the flux. For this source, the two model components listed are simply for each galaxy in the pair, rather than a 
multicomponent fit to a single galaxy. The fit to the compact galaxy was held constant with a de Vaucouleurs profile, while models were allowed to vary for the 
second, extended galaxy. Something close to an exponential or low Sersic index gives the best fit for this extended galaxy. A de Vaucouleurs profile appears to 
be ruled out by the ACS data. 



Inspection of the HST images for B2 1018+34 ER2 reveals 
a striking grand design spiral arm structure. The arms are 
clearly bluer than the central region, especially a few notice- 
able knots of star formation in the ACS image. Some of 
these knots are also visible in the NICMOS image, notably 
the knots that appear to be on the ends of a possible bar, or the 
innermost region of the spiral arms. In Paper I we concluded 
that this galaxy consisted of a younger stellar population than 
the rest of the galaxies in our sample, and noted the presence 
of [O II] emission in the rest-frame UV spectrum. Asymme- 
tries such as spiral arms are not accounted for in the GALFIT 
modeling, however, a smooth exponential profile with the ad- 
dition of a central bulge that contributes 24% of the total light 
fits the underlying structure of the galaxy well (Pig. [5j. After 
subtraction of the model galaxy, arm structure that was not as 
evident in the NICMOS image as it was in the ACS image 



becomes clearly visible in the residuals. 

4.4. TXS 121 1+334 ER1 

TXS 1211+334 ER1 appears to be a somewhat elongated 
elliptical galaxy. The best fit GALFIT profile is a Sersic 
n = 3.90 with a small effective radius (Fig. 01. An expo- 
nential profile with a small effective radius plus a central PSF 
also fits the central regions of the galaxy, but is too low at 
large radii. It is possible that this galaxy is still in the pro- 
cess of relaxation into a classical elliptical galaxy. The resid- 
uals from the Sersic fit show there is excess flux along the 
axis of elongation (see FigQ]and residual images in Fig. [6j>, 
possibly remnants of a much earlier merging event. If there 
was a major merging event in the past, it could not have pro- 
duced much star formation since the age of the dominant stel- 
lar population, as shown by both broad-band photometry and 
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FIG. 4. — Surface brightness profile for J094258+4659 ERO. The best fit model is an exponential disk with an r 1 / 4 bulge that contributes 32% of the total light 
at rest-frame R-band and 25% of the light at rest-frame U-band. 



medium-resolution spectra (Paper I), is extremely old (> 1.4 
Gyr), with little evidence for any admixture from a younger, 
UV-dominated population. The ACS profile is in agreement 
with the NICMOS profile, also favoring a classical de Vau- 
couleurs elliptical galaxy. 

4.5. 4C 15.55 ER2 

From the HST imaging, 4C 15.55 ER2 appears from sim- 
ple visual inspection to be either an elongated elliptical or an 
inclined disk galaxy. Previous results from high S/N ground- 
based AO imagi ng concluded that th e galaxy was a nearly 
pure disk galaxy dStockton et alj200q) with a small point-like 
core. The NICMOS imaging corroborates this result. The best 
GALFIT model is an exponential profile with the addition of a 
central PSF that contributes 5% of the total light (Fig. 13). A 
Sersic profile with an index, n = 1.5 also fits the data well. 
In either case, it is clear that this galaxy is dominated by a 
massive disk of old stars. The ACS profile yields a similar 
result, with a somewhat smaller Sersic index (n = 1.35) or 
an exponential disk with a slightly smaller core (contributing 



2.5% of the total light). An r 1 / 4 -law profile clearly provides a 
much worse fit to both the ACS and NICMOS imaging data. 

4.6. TXS 1812+412 ER2 

TXS 1812+412 ER2 appears to be fairly round with some 
extended diffuse emission around the core. This fainter emis- 
sion is nearly symmetrical in both the ACS and NICMOS im- 
ages. The best fit GALFIT profile is a combination of some- 
thing like a disk (an exponential profile or a Sersic profile 
with a low index) plus a small core (either point-like, or a 
small de Vaucouleurs bulge). Figure [8] shows the fits from 
different single and multi-component models. Pure exponen- 
tial and de Vaucouleurs profiles do not provide a good fit to 
the data as can be seen most clearly in the difference panels at 
the bottom of Figure [8] 

5. COLOR GRADIENTS 

Elliptical galaxies in the local universe exhibit color gra- 
dients, with t he inner regions being redder than the outer- 
most regions dFranx et al.1 fl989t iPeletier et al.1 [1990). This 
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FIG. 5. — Surface brightness profile for B2 1018+34 ER2. The best fit model is an exponential disk with an r 1 / 4 bulge that contributes 24% of the total light 
at rest-frame R-band and 12% of the light at rest-frame U-band. Spiral arms contribute to radial asymmetries that are difficult to model and are most readily 
apparent in the ACS profile. 



trend has been shown to be largel y due to a metallicity gra- 
dient, r ather than an age gradi e nt (jDavies. Sadler. & Peletierl 
1993b iKobavashi & Arimotol Il999t iTamura et all 120001: 
Hinklev & Iml 12001 ). The strong gravitational potential in 
the center of a galaxy contributes to the retention of high- 
metallicity gas, whereas in the outer regions of a galaxy, 
stellar winds from supernovae are efficient in expelling this 
metal-rich gas, leading to the observed color gradient. Such 
a gradient can be explained quite naturally with a monolithic 
collapse formation scenario, in which the bulk of the stars are 
formed in a single short burst and passively evolve thereafter. 
Metallicity gradients in this scenario are likely in place shortly 
after formation, and remain roughly unchanged as the galaxy 
evolves. On the other hand, hierarchical formation scenar- 
ios predict that massive ellipticals form over longer periods of 
time, gradually being built up from smaller sub-units. These 
models predict that elliptical galaxies consist of several dis- 
tinct stellar populations, with mergers driving gas and dust 



into the central regions promoting mixing and fueli ng star 
forma tion which in turn dilutes color gradients (e.g., IWhitd 
1980), possibly even leading to inverted color gradients (e.g., 
Im et alJuOOll) . At high z, color gradients can be used to infer 
formation mechanisms since the galaxies are closer to their 
formation epoch and the youngest stars will be significantly 
bluer than the oldest stars. 

We have examined color gradients in our seven EROs to 
better understand what possible age or metallicity gradients 
are present. In order to do this we first needed to resample 
the NICMOS and ACS images to put them on the same pixel 
scale, and then convolve them with the PSF from the other 
instrument. We used the same PSFs that were generated for 
the morphological modeling with the exception that the PSFs 
were resampled to the proper pixel scale before convolution. 
This process eliminates any artificial gradient due to PSF dif- 
ferences between the F8 14W and F160W images. We then ro- 
tated the images using header keywords that defined the tele- 
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FIG. 6. — Surface brightness profile for TXS 121 1+334 ER1. The best fit model is a Sersic profile with an index n - 
The two models overlap at nearly all radii. 



3.90, or a classic de Vaucouleurs profile. 



scope roll angle for each observation, and shifted them using 
the galaxy centroids as a primary alignment. When possible, 
other stars in the field were used to double check the rota- 
tion and offsets in order to ensure that differences in galaxy 
centroids between the blue and red components did not affect 
the alignment. The background was subtracted from both im- 
ages before producing the final color maps. We have made 
postage -stamp cut-out images of the galaxies using the SEx- 
tractor (iBertin & Arn outs 1996) segmentation maps in order 
to mask out the surrounding noisy background, which fluctu- 
ates around zero. These segmentation maps define the galaxy 
regions based on all contiguous pixels above a flux threshold 
of 2.5 a in the NICMOS imaging. In terms of effective radii, 
r e , the color maps generally include all galaxy flux within 
~2r e , except for TXS 0145+386 ER1, ER2, and B2 1018+34 
ER2 where this flux threshold retained all galaxy flux within 
1.0 to 1.5 r e . These color maps are shown in Figure|9] 

The two galaxies that clearly follow r 1 / 4 -law profiles, TXS 
0145+386 ER1 and TXS 1211+334 ER1, have relatively flat 



colors, consistent with little or no color gradients. These 
galaxies are uniformly red at least to ~ 1 — 2r e , where the 
S/N becomes too low to precisely determine any further color 
variation. There is no evidence for any inverted color gradi- 
ent, such as th at seen in roughly 50% of the Hubble Deep 
Field galaxies (Me nanteau et alJ l2001h . which implies that 
these galaxies most likely have not suffered any recent star 
forming events. This conclusion is corroborated by the in- 
dividual ACS and NICMOS images, which show extremely 
relaxed systems with similar morphologies both short-ward 
and long-ward of the 4000A-break. However, the flatness of 
the color gradients in these two elliptical galaxies is quite dif- 
ferent from what we observe in massive ellipticals in the local 
universe. It is possible that some amount of dry merging has 
already occurred in these galaxies, diluting any initial metal- 
licity gradients present from formation. It is interesting to 
note that a "companion" galaxy to the east of TXS 0145+386 
ER1 has the same color as that of ER1 . This companion could 
be indicative of a possible future minor merger event. 
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FIG. 7. — Surface brightness profile for 4C 15.55 ER2. The best fit model is an exponential disk with a small point-like core that contributes between 2.5 — 5' 
of the total light. A Sersic profile with a low index (n = 1.49) is also a good fit. A pure r 1 / 4 -law profile does not provide a good fit to either the NICMOS i 
ACS observations. 



The color map of TXS 0145+386 ER2, which consists of 
another merging pair, shows that one galaxy is clearly much 
redder than the other. The observed color gradient in the main, 
red galaxy is most likely too steep to be reproduced by metal- 
licity alone. Furthermore, a blue tail appears to fan out to the 
northwest from the compact source, which could indicate a 
slightly younger disk component disrupted by the merger. The 
combined colors of this merging pair yield an SED that is con- 
sistent with extremely little or no dust (Paper I). If one of the 
two galaxies has a significantly younger population of stars, 
it would imply a much older age for the dominant stellar pop- 
ulation in order to match the combined SED. Rest-frame UV 
spectroscopy, which was not obtained for this source, would 
be extremely helpful in determining whether any resulting star 
formation, however small, is occurring due to this merger, or 
whether it truly is a "dry" merger. 

B2 1018+34 ER2, which was determined in Paper I to be 
a younger, possibly star-forming galaxy has a very interest- 
ing color profile. The central bulge is red, and a dust-lane is 
visible running in front of the bulge. Bluer regions from the 



spiral arms visible in the ACS image are apparent on either 
side of the bulge. Overall, it seems clear that this is a clas- 
sic grand design spiral galaxy, with star formation occurring 
more recently in the spiral arms than in the central bulge re- 
gion. However, the possible presence of [O III] contaminating 
the J-band photometry (Paper I), may indicate that an AGN 
is partially responsible for the observed color gradient. 

J094258+4659.2 ERO, discovered bv llve et alj d2000b . also 
shows an interesting color profile. Results from the surface 
brightness profile fitting for this galaxy concluded that it was 
dominated by an exponential disk, with a bulge component 
that contributed ^30% of the light. Here we can see clearly 
that one side of the galaxy is bluer than the other, probably 
indicating the presence of dust on the southeast side. The 
bulge is clearly red, and its stars are most likely slightly older 
than those in the disk. If we take the color difference between 
the NW and SE sides of the disk and assume the difference 
arises solely from reddening effects, the extinction we infer 
based on the standard Galactic extinction curve (Ry = 3) is 
consistent with the globally averaged extinction derived from 
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FIG. 8. — Surface brightness profile for TXS 1812+412 ER2. The best fit for this galaxy is a two component model of either a Sersic profile with a low index 
(n = 1.50) plus a point-like core that contributes 13% of the light, or an exponential profile plus a small core that contributes ~16% of the light. Disks with 
compact r 1 / 4 -law bulges (not shown) provide nearly as good a fit as the stellar cores. Pure exponential and pure de Vaucouleurs models are ruled out. 

star formation occurring in this galaxy. All signs indicate that 
this is an extremely relaxed, massive disk of old stars that 
must somehow evolve into something resembling a massive 
elliptical galaxy by the present day. 



the SED fitting jive et al.ll2003l) . 

TXS 1812+412 ER2 was determined to have a best-fit 
Sersic profile with n = 1.5 plus either a small point-like or 
compact r 1 / 4 -law core. It appears to have a rather classic 
elliptical color gradient, however, with a red central region, 
becoming bluer in the outermost regions. This mild gradi- 
ent could also be explained by a disk component that is only 
mildly younger than the central bulge. 

Lastly, and perhaps most interesting, is the color profile of 
4C 15.55 ER2. This galaxy is quite clearly dominated by an 
old stellar population (1.8 Gyr; Paper I) and is also clearly 
dominated by an exponential profile (Fig. |7J. The observed 
color gradient in this galaxy, however, is reminiscent of that 
found in local giant elliptical galaxies. This gradient is al- 
most certainly due to a metallicity gradient, rather than an age 
gradient, given the fact that something very close to a single 
instantaneous burst of star formation appears to be required 
in order to fit the photometry. There is little evidence from 
either the photometry or the spectroscopic data for any recent 



6. DISCUSSION 

6.1. Massive Disks of Old Stars 

Morphological classes for EROs at high z have largely been 
assumed to be correlated with their spectral type. That is, 
dusty starbursts are generally assumed to be disk-dominated, 
while galaxies with evolved stellar populations are assumed 
to be bulge-dominated. We have found three cases, however, 
or roughly half of our sample, that exhibit exponential (or 
nearly exponential) profiles and are dominated by extremely 
old stars, plus a fourth ga laxy with an exp onential profile and 
somewhat younger stars. lYan et all (120041) find a similar result 
for their sample of EROs at z ~ 1, with a little over half (7/12) 
of their absorption-line-only syst ems being disk-dominated 
(see also. lYan & Tho mpson 2003). In fact, even among their 
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TABLE 3 
Properties of Old Galaxies 



Galaxy 


2 


best fit Sersic model 


two-component model 


r eff 


class 










(kpc) 




TXS 0145+386 ER1 


1.4533 


n = 4.3 




6.65 


spheroid 


ER2 


1.459* 


n = 1.9 &n = 4.0 




5.12 & 0.77 


merger 


J094258+4659.2 ERO Rl 


1.50* 


n = 3.6 


expdisk + 32% bulge 


7.34 


disk 


B2 1018+34 ER2 


1.4057 


n = 3.6 


expdisk + 24% bulge 


5.95 


disk 


TXS 121 1+334 ER1 


1.598 


n = 3.9 




2.22 


spheroid 


4C 15.55 ER2 


1.412 


n = 1.5 


expdisk + 5% PSF 


3.15 


disk 


TXS 1812+412 ER2 


1.290 


n = 1.5 + 13% PSF 


expdisk + 16% PSF 


3.63 


disk 



* The best-fit photometric redshift is given as no spectroscopic redshift is available for these two sources. 



galaxies that show small amounts of [O II] emission, they 
note that the underlying stella r populations often ap pear to be 
quite old. On the other hand, Abr aham etail d2007l) find that 
nearly all of their spectroscopically quiescent systems in the 
Gemini Deep Deep Survey are bulges, and on the basis of this 
conclude that morphology is a reasonable indicator of spectral 
activity. Their approach to morphol ogical classification, how- 
ever, is quite different from that of lYan & Thompson] (120031) 
or from our own approach. They use a technique based on an 
asymmetry index and a Gini index, which can be though t of as 
a generalized concentration index (Abraham et al. 2003]). The 
requirement for a galaxy to be classified as bulge-dominated 
under this system is a low asymmetry (< 0.15) and a high 
Gini index (> 0.50). This assumes that disk galaxies will be 
more asymmetric due to spiral arms and knots of star forma- 
tion, and less centrally concentrated than their elliptical coun- 
terparts. However, a disk of old stars can be quite symmetric, 
so there is reason to be cautious when using this classifica- 
tion system. Furthermore, it is well known, as Abraham et al. 
point out, that the Asymmetry-Gini (AG) classification cannot 
distinguish between E and SO galaxies. 

In fact, we have found that the inability of the AG method 
to distinguish between morphological classes appears to ex- 
tend to galaxies of all types that are dominated by old stars. 
J094258+46592 ERO, 4C 15.55 ER2, and TXS 1812+412 
ER2, which appear to be dominated by disks of old stars from 
our NICMOS surface brightness profile fitting, all have asym- 
metry indices near 0.05 and Gini indices near 0.50 (two of 
them having a Gini index slightly greater than 0.50 and one 
being slightly l ess than 0.50), as de termined by the MOR- 
PHEUS 7 code dAbraham et all 120071) . These values would 
place them into the lAbraham et al.l "early type" category, 
which is clearly a misclassification in terms of morphology. 
We therefore have based our morphologies on radial surface 
brightness profiles, since this method does not assume any- 
thing about the distribution of starlight in the galaxies we are 
measuring, and appears to give the most robust results, espe- 
cially when dealing with evolved stellar populations. 

For the four galaxies in our sample that were determined 
to have strong exponential components to their radial light 
profiles, we examined their sizes and compare them to those 
of present-day disk galaxies. J094258+4659 ERO and B2 
1018+34 ER2 have effective radii of 7.3 and 6.0 kpc, respec- 
tively, which correspond to disk scale lengths of 4.4 and 3.6 
kpc. This is roughly consistent with average disks in the lo- 
cal universe (Ivan der Kruitlll987h which have scale lengths of 

1 MORPHEUS is available at http://odysseus.astro.utoionto.ca/~abraham/Moipheus/ 



5. 5 ±2 .7 kpc, and i s comparable to the massive disks found by 
lLabbe et all (120031) at z > 2. These galaxies could well be the 
passively evolved descendants of such a population. The ef- 
fective radius for 4C 15.55 ER2, on the other hand, is only 3.2 
kpc, corresponding to a scale length of 1 .9 kpc, and the effec- 
tive radius of TXS 1812+412 ER2 is 3.6 kpc, corresponding 
to a scale length of 2.1 kpc. While the number statistics are 
low, taken together there is a general trend that these old disk 
galaxies at z ~1.5 are 45% smaller than present-day disks. 
This is in line with recent work on t he size evolution predi cted 
by CDM dark matter halo growth dSomerville et alj|2008l) . 

While their size may be well explained under stan- 
dard CDM predictions, the three oldest disk galaxies, 
J094258+46592, 4C 15.55 ER2, and TXS 1812+412 ER2, 
pose an interesting problem for hierarchical formation sce- 
narios nonetheless. 4C 15.55 ER2, in particular, is difficult to 
explain since it appears to be a nearly pure disk of 2 x 10 11 M Q 
that completed star formation by z > 2.4. Any mergers that 
may have occurred between 1.4 < z < 2.4 would have to 
have been minor mergers of gas poor systems in order to pre- 
serve both the disk structure and the age of the stellar popu- 
lation. It therefore seems unlikely that mergers could be the 
dominant formation mechanism for such a galaxy. Instead, 
something more similar to monolithic collapse, where star- 
formation proceeds rapidly and efficiently in a single burst 
spread more-or-less evenly across the disk, seems much more 
reasonable. 

The selection of our sources near radio-loud quasars likely 
places them in more overdense regions that will eventually 
form present-day clusters. As such, they are unlikely to be 
representative of all quiescent galaxies at z ~ 1.5. How- 
ever, the galaxies were chosen purely by photometric criteria, 
and there is no reason to believe that disks would be pref- 
erentially selected. An equivalent selection in the local uni- 
verse would result almost exclusively in elliptical and SO (i.e., 
bulge-dominated) galaxies. Therefore, while not necessarily 
representative of all quiescent galaxies at high redshift, it is 
clear that at least some massive galaxies are unlikely to have 
assembled the bulk of their stars through mergers. 

6.2. Size Evolution 

A number of recent studies have found evide nce for highly 
compact, massive galax ies at z > 1.4 (e.g.. iTruiillo et al. 
2006; ZirmetalJ l2007t iToft et all l2007t IvanDokkum et al. 
2008). However, the majority of our sources are not espe- 
cially compact when compared with local galaxies of similar 
mass. The typical effective radius of a 2 x 10 11 M W galax y 
in the SDSS is 6 kpc at rest-frame I-band dShen et al.ll2003l) . 
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FIG. 9. — Color maps for seven old galaxies from HST imaging. TXS 0145+386 ER2 consists of two galaxies; the centroid for the smaller galaxy in the 
merging pair is marked by a black cross. Scale for each image is denoted by a white bar which represents 1" in angular size. Color wedges on the right hand side 
of each panel represent the F814W— F160W color scale for each corresponding galaxy in units of magnitude per pixel. 

The galaxies presented here are typically ^6 kpc at rest-frame 
R-band. Three galaxies (TXS 121 1+334 ER1, 4C 15.55 ER2, 
and TXS 1812+412 ER2) are more compac t than this by a 
factor of ~2. Note, however, that we use toe IChabriei] 42003) 
IMF which results in mass estimates that are a factor of 1.5 
times smaller than mas ses calculated using the lKroupald2001l) 
IMF as in the SDSS (lLonghetti et al.1 120051) . If we apply 
this correction factor to our mass estimates, then the typical 
size for galaxies of similar mass in the SDSS would be 7.0 
kpc, and our galaxies would generally be 15 % smaller than 
this, with three more than 50% smaller than this. There is a 
slight tendency toward smaller effective radii at longer wave- 
lengths (e.g., our rest-frame i?-band vs. rest-frame [/-band 
data), which could increase the discrepancy further if we had 
morphological information at rest-frame /-band to compare 
with the SDSS. However, in most cases size appears to be 
roughly independent of observed wavelength, as evidenced 
by the approximately flat color gradients seen in many of 
these galaxies. None of our galaxies are as compact a s the 
r < lkpc galaxies found at z > 2 by iToft et al.1 (120071) and 
Ivan Dokkum et a l. (2008). Our three s mallest galax i es are 
consistent with galaxies at z ~ 1.4 in the lTruiillo et al.l d2006l) 
sample, while the rest are significantly larger. 



It is interesting to note that the galaxies with the lowest 
dust content (Ay < 0.6) from the photometric SED fitting 
(with the exception of TXS 0145+386 ER2, which appears 
to be undergoing a dry merger), also appear to be the most 
compact, regardless of the disk or bulge-like nature of their 
morphology. This may be a further manifestation of the rela- 
tionship between star- formation and siz e at z ~ 2.5 found by 
IZirm et alj d2007l) and IToft et all d2007l) . All of our galaxies 
are best fit by single, instantaneous bursts of star formation 
at high redshift rather than constant star formation models 
an d would therefore all fall under the " quiescent" category 
of IZirm et all d2007l) and lToft et ail d2007l) . However, the evi- 
dence for dust in galaxies with larger effective radii (and like- 
wise, the lack of dust in the most compact sources) may point 
to a relationship between star-formation history and size. In 
this scenario galaxies that form rapidly and undergo intense, 
highly efficient star formation early-on end up as compact ob- 
jects, either disks or bulges, while galaxies where star forma- 
tion is slightly less efficient may leave behind more dust and 
may form more gradually, leading to extended morphologies. 
While the number statistics are low, this idea would appear 
to be supported by the fact that galaxies wit h similarly aged 
stellar populations at higher redshifts (e.g., Toft et al. 2007; 
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IStockton et all 120081: Ivan Dokkum et all 120081) tend to have 
even more compact morphologies. These z > 2 galaxies must 
have formed over even shorter time periods than the galaxies 
in the present z ~ 1.5 sample. 

6.3. Role of "Dry Mergers " 

Regardless of how these galaxies formed originally, it is 
clear that given the age of the universe, over the remaining 
~12 billion years from the end of star formation until the 
present day, it is likely that they will undergo a number of 
subsequent mergers. This is even more evident for the disks 
of old stars that we see at z ~ 1.5, since we do not find any 
comparable galaxies in the local universe. Some mechanism 
must be responsible for converting them into objects that more 
closely resemble either SO or giant elliptical galaxies in the 
present-day. Furthermore, the lack of luminous blue galaxies 
in the local universe requires that this mechanism must not 
destroy or dilute the old stellar populations present in these 
galaxies. 

There has been much discussion recently about "dry merg- 
ers", in which little or no cold gas is available to form new 
stars, as an important mechanism for preserving old stel- 
lar populations and ge nerating the mass i ve spheroids we see 
in the local univ e rse feell et alj [2004t Ivan Dokkum] 120051: 
iFaber et alj[2007l) . iBell et alj d2004l) find that the rest-frame 
B-band luminosity density of red galaxies evolves only mildly 
between < z < 1.1. This requires at least a factor of two 
increase by z = in the stellar mass density present on the 
red sequence in order to compensate for dimming d ue to con- 
tinued passive evolution of the stellar populations. IBell et all 
attribute this evolution to both gas poor mergers and to trunca- 
tion of star formation in blue star-forming galaxies. IBell et all 
(2006) also studied the dry merger fraction among galaxies in 
the GEMS survey and found that luminous early-type galax- 
ies on the red sequence will undergo on average ~1 major dry 
merge r between z < 0.7 and the present day. van Dokkum 
(2005) finds a similar result, with ^70% of bulge-dominated 
galaxies experiencing a dry merger with a median mass ra- 
tio 1 : 4 in the recent past (z < 0.2). He derives a current 
mass accretion rate of roughly 10% Gyr" 1 , which again im- 
plies about a factor of 2 increase in stellar mass density in 
luminous red galaxies between < z < 1. 

TXS 0145+386 ER2 may be an example of this dry merging 
process in action. The colors of the galaxies in this merg- 
ing pair are consistent with that of old stellar populations 
rather than dusty starbursts (Paper I), and there is no evidence 
for any significant amount of recent star formation. Recall 
that the compact, northwest galaxy in this pair was best fit 
by a very small effective radius, while the second, more ex- 
tended galaxy to the southeast appeared to be best fit by some- 
thing similar to an exponential disk, especially at rest-frame 
U-band. In addition, the color map for this galaxy showed 
evidence for some possibly younger (although still relatively 
old) stars fanning out from the compact NW galaxy, which 
may imply that this spheroid is the central bulge of a disk 
galaxy being torn apart by the interaction. If these scenar- 
ios are correct, this could imply a dry merger between two 
old disk galaxies, which would provide the much needed link 
between these oddities at z = 1.5 and present-day massive 
ellipticals. 

Another of our galaxies, TXS 121 1+334 ER1, was best fit 
by an r 1 / 4 -law, but residuals along the axis of elongation sug- 
gest a more complex formation history. It is possible that this 
galaxy suffered a merger in its recent past. However, given 



that no signatures of star formation are visible in the rest- 
frame UV spectrum of the source and that the SED favors 
a stellar population that is at least 1 .4 Gyr old (Paper I) with 
a current star formation rate of < 0.2 M Q yr _1 (Table [2]), 
any merging event must have been gas poor. Taken together, 
TXS 0145+386 ER2 and TXS 1211+334 ER1 may represent 
different stages of a similar process. 

7. SUMMARY AND CONCLUSIONS 

We have studied a sample of seven evolved galaxies at 
z ~ 1.5 with high resolution imaging from HST These are 
summarized below. 

1. TXS 0145+386 ER1 is a relaxed elliptical galaxy at 
z = 1.4533 with a dominant stellar population that 
is ~1 Gyr old (Paper I). There is no evidence for any 
color gradient across the galaxy, implying the stars are 
uniformly old, yet may have been mixed through dry 
merging events in the past. A "companion" galaxy 2" 
to the east has the same color as ER1 and may indicate 
a current or future minor merging event. There are no 
visible signs of interaction between these two galaxies, 
but these features may be too faint to observe at this 
redshift, especially if they are dominated by old stars. 

2. TXS 0145+386 ER2 is a pair of merging galaxies at 
z ~ 1.46. One galaxy is extremely compact, while 
the other is much more extended and redder in color. 
Their combined flux gives an SED of an extremely old 
stellar population (2.5 Gyr) with no dust (Paper I). This 
result implies that no recent star formation is occurring 
and that they are undergoing a dry merger. The ACS 
image of this pair shows some irregular structure in the 
extended galaxy which, if due to small amounts of star 
formation and/ or dust, implies even older ages for the 
majority of the stars in order to produce an average age 
of 2.5 Gyr. There is some evidence from the color map, 
as well as a visible elongation in the ACS image, for 
bluer stars in a tail outward from the compact galaxy to 
the northwest; however the S/N is low in this outermost 
region. 

3. J094258+4659.2 ERO is a disk-dominated galaxy at 
z ~ 1.5 with a bulge that contributes roughly 30% 
of the total light. Th e dominant stel lar population is 
on average 1 Gyr old dive et a l. 2003). A color gradi- 
ent is present, with the central bulge being redder than 
the disk, and the disk is noticeably bluer on the north- 
west side and redder on the southeast side. This is most 
likely due to slightly different stellar populations in the 
bulge and disk, as well as dust obscuration along the 
southeast side of the disk. 

4. B2 1018+34 ER2 is a grand design spiral galaxy at z = 
1.4057 with a stellar population that is on average 0.64 
Gyr old (Paper I), although it is likely that the majority 
of the stars are much older (~ 1.9 Gyr, Table|2} with a 
small contribution by mass from younger stars. Spiral 
arms are most clearly visible in the ACS image, while 
the NICMOS image has a smoother profile. A color 
map of the galaxy shows that a dust lane runs from the 
southeast to the northwest in front of the bulge. 

5. TXS 121 1 +334 ER1 appears to be a classic elliptical 
galaxy at z = 1.598 with an age of 1.4 Gyr (Paper I). 
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Residuals from subtracting an r^-law profile from the 
NICMOS image indicate that the galaxy may have suf- 
fered a merging event in the past and may still be in the 
process of relaxing. There is no evidence for any sig- 
nificant color gradient across the galaxy, which may be 
due to population mixing from prior merging events. 

6. 4C 15.55 ER2 is essentially a pure disk of old (1.8 Gyr, 
Paper I) stars at z — 1.412. A moderate color gradient 
is visible in this galaxy, with the central regions slightly 
redder than the outermost regions, which could indicate 
a metallicity gradient. It is important to note that there 
does not appear to be a significant bulge component to 
this galaxy (< 5% of the total light), so this color gra- 
dient cannot be explained by a difference in stellar pop- 
ulations between a bulge and disk component. 

7. 7X5 1812+412 ER2 appears to be another disk- 
dominated galaxy at z = 1.290 with an age of 1.4 Gyr 
(Paper I). The best-fit profile for this source requires a 
~ 10 — 15% contribution from a compact core. There 
is a mild color gradient visible for this source, with the 
central regions being slightly redder than the outer re- 
gions. 

From the morphologies of these galaxies it is clear that 
age, or spectroscopic activity is not necessarily a surrogate for 
morphological class, and vice versa. Three of these galaxies 
(J094258+4659.2, 4C15.55 ER2, and TXS 1812+412 ER2) 
have strong exponential components to their radial light dis- 
tributions, and both their broad-band SEDs and r est-frame 
UV c ontinua point to evolved stellar populations (live et alJ 
2003; Paper I). The existence of massive evolved disks at 
high redshift is a challenge for current semi-analytic models 
of galaxy formation, since it implies that they have not under- 



gone a single major merging event since their formation more 
than 1 Gyr earlier. Instead, they appear to have formed by 
a mechanism more similar to monolithic collapse. Further- 
more, their existence is problematic for new theories that rely 
on AGN feedback in galaxy formation in order to reproduce 
the observed large scale properties of galaxies, since power- 
ful AGNs are likely triggered by major mergers and prefer- 
entially exist in luminous spheroids. It would seem that an- 
other mechanism must be responsible either for triggering the 
AGN feedback or for terminating star formation in these disk 
galaxies at early times, possibly by the rapid formation and 
efficient processing of all available gas into stars in a single 
burst of activity. 

Dry mergers may be the important link between these old 
galaxies and present-day massive ellipticals. We have found 
at least one example of this process in action. Dry mergers ap- 
pear to be an important evolutionary phase in most luminous 
red galaxi es at z < 1 dBell et alJl200H Ivan Dokkuml 120051: 
iBell et alJl2006h . However, massive disks of old stars that 
have not suffered a single major merger are present at early 
times, implying that multiple formation mechanisms may be 
required to explain how the earliest massive galaxies formed. 
It will be important to continue mo rphological studies of old 
galaxies at redshifts z > 1.5 (e.g JStockton et aTl l2008i) in or- 
der to determine how common massive disks of old stars are 
in the early universe. 
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